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Abstract. Many contemporary applications, for example, cataloging of galax- 
ies, document analysis, face recognition, learning theory, image processing, op- 
erate with a large amount of data which is often represented as a graph embed- 
ded into a high dimensional Euclidean space. The variety of problems arising 
in contemporary data processing requires development on graphs such topics 
of the classical harmonic analysis as Shannon sampling, splines, wavelets, cu- 
bature formulas. The goal of the paper is to establish cubature formulas on 
finite combinatorial graphs. The results have direct applications to problems 
that arise in connection with data filtering, data denoising and data dimension 
reduction. 
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1. Introduction 

Many contemporary applications, for example, cataloging of galaxies, document 
analysis, face recognition, learning theory, image processing, (see [8] for references) 
operate with a large amount of data which is often represented as a combinatorial 
graph embedded into a high dimensional Euclidean space. In such situation any 
processing of data depends of ones ability to perform Harmonic analysis of functions 
(signals) defined on graphs and not on the regular Euclidean spaces. By Harmonic 
analysis on graphs we understand not only the traditional analog of Fourier analysis 
as the spectral analysis of the Laplace-Beltrami operator. The variety of problems 
arising in contemporary data processing requires development on graphs such clas- 
sical themes as Shannon sampling, splines, wavelets, cubature formulas. All these 
topics are definitely underdeveloped in the setting of graphs. The goal of the paper 
is to establish cubature formulas on finite graphs. The present article utilizes our 
previous work on Shannon sampling and splines on combinatorial graphs [3]- [7]. 
Our results can find applications to problems that arise in connection with data 
filtering, data denoising and data dimension reduction. 

Given a real- valued function / on a set of vertices V of a combinatorial graph its 
"integral" is J2vev f( v )- O ur § oa ^ ^ s to develop a set of rules (cubature formulas) 
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which allow for approximate or exact evaluation of integrals by using values of a 
function on a subset U C V of vertices. 

In section 2 we review our results [1] about variational interpolating spines on 
graphs and describe an algorithm which allows an effective computation of vari- 
ational splines. In section 3 by using interpolating splines we develop a set of 
cubature formulas. Theorem 3.1 shows that these formulas are exact on the set of 
variational splines. Theorem 13.21 provides an approximate cubature formula for all 
functions on the graph and Theorem l3.4l shows that this cubature formula is asymp- 
totically exact for banlimited functions of low bandwidths. At the end of section 
3 it is proved that cubature formulas based on interpolating splines are essentially 
optimal (Theorem I3.5j) . 

As it is explained in the Example after Theorem 13.41 for a cycle graph of 1000 
vertices a set of about 670 "uniformly" distributed vertices is sufficient to have 
asymptotically exact cubature formulas for linear combinations of the first 290 
eigenfunctions (out of 1000) of the corresponding combinatorial Laplace operator. 

It is worth to note that all results of section 3 which provide errors of approxi- 
mation of integrals of functions on V through their values on a U C V reflect 

1) geometry of U which is inherited into the quantity \/\V\ — \U\ = y/\S\ and 
into the Poincare constant A (see section 3 for definitions), 

2) smoothness of functions which is measured in terms of combinatorial Laplace 
operator. 

In section 4 we develop a different set of cubature formulas which are exact on 
appropriate sets of bandlimited functions. The results in this section are formulated 
in the language of frames and only useful if it is possible to calculate dual frames 
explicitly. Since in general it is not easy to compute a dual frame we finish this 
section by explaining another approximate cubature formula which is based on the 
so-called frame algorithm. As the formula (|4.24[) shows the rate of convergence of 
this cubature rule depends solely of the way subset U (from which all the samples 
are taken) is embedded into entire graph V. To be more specific it depends on the 
ratio of the number K (U) of "in going" edges to the number D (U) of "out going" 
edges. 



2. Variational (polyharmonic) splines on graphs 

We consider finite or infinite and in this case countable connected graphs G = 
(V(G),E(G)), where V(G) = V is its set of vertices and E(G) = E is its set of 
edges. We consider only simple (no loops, no multiple edges) undirected unweighed 
graphs. A number of vertices adjacent to a vertex v is called the degree of v and 
denoted by d(v). We assume that degrees of all vertices are bounded from above 
and we use notation 

d(G) = ma,xd(v). 

The space L2(G) is the Hilbert space of all real valued functions / : V — > K with 
the following inner product 



(/,<?> = £ f{v)g(v) 

vEV 
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and the following norm 

11/11 = ll/llo= (ei/wi 2 

\vev 

Let A be the adjacency matrix of G and I? be a diagonal matrix whose entree 
on main diagonal are degrees of the corresponding vertices. Then we consider the 
following version of the discrete Laplace operator on G 

(2.1) C = D-A, 

or explicitly 

Cf(v) = J2(f(v)-f(u)),feL 2 (G), 

where notation u ~ v means that u and v are adjacent vertices. 

Operator(matrix) C is symmetric and positive definite. Let E w (£) be the span of 
eigenvectors of C whose corresponding eigenvalues are < uj. The invariant subspace 
E w (£) is the space of all vectors in L2(G) on which £ has norm w. In other words 
/ belongs to E w (C) if and only if the following inequality holds 

(2-2) \\C s f\\<w s \\f\\,s>Q. 

For a given set of indices / (finite or infinite) the notation l 2 will be used for the 
Hilbert space of all sequences of real numbers y — {y u }, for which |y„| 2 < oo. 

Variational Problem 

Given a subset of vertices U = {u} C V, a sequence of real numbers y = {y u } S 
l 2 ,u G U, a natural k, and a positive e > we consider the following variational 
problem: 

Find a function Y from the space L 2 (G) which has the following properties: 

1) Y{u)=y u ,ueU, 

2) Y minimizes functional Y — > ||(e/ + £) fe F||. 

We show that the above variational problem has a unique solution Y^ v . 

For the sake of simplicity we will also use notation Y£ assuming that U and e 
are fixed. 

We say that Y% is a variational spline of order k. It is also shown that every 
spline is a linear combination of fundamental solutions of the operator (el + C) k 
and in this sense it is a polyharmonic function with singularities. Namely it is 
shown that every spline satisfies the following equation 

(2.3) [eI + C) 2k Yy = J2^u6 u , 

ueu 

where {a u } ue u — < a u (Y^) > is a sequence from l 2 and S u is the Dirac measure 

at a vertex u e U. The set of all such splines for a fixed U C V and fixed 
k > 0, e > 0, will be denoted as y(U 7 k, e). 

A fundamental solution F£ (= F% ,u e V, of the operator (el + C) k is the 

solution of the equation 




(2.4) 



(eI + C) k F^ = S u , fceN, 
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where S u is the Dirac measure at u £ V(G). It follows from ()2.3|) that the following 
representation holds 

= £a u i^. 

v.eu 

It is shown in [4] that for every set of vertices U = {u}, every natural k, every 
e > 0, and for any given sequence y — {y u } £ h, the solution of the Variational 
Problem has a representation 

(2-5) l?=£w» L *> 

where L k is the so called Lagrangian spline, i.e. it is a solution of the same Varia- 
tional Problem with constraints L k (v) = 5 u _ Vl u £ U, where 5 U>V is the Kronecker 
delta. It implies in particular, that y(U,k,e) is a linear set. 

Given a function / 6 -La(G0 we will say that the spline yI interpolates / on U 
if y/(u) = f( u ) for all u £ U. 

Algorithm for computing variational splines. 

The above results give a constructive way for computing variational splines. 
Suppose we are going to construct splines which have prescribed values on a subset 
of vertices U C V. 

1. One has to solve the following \U\ systems of linear equations of the size 

\V\ x \V\ 

(2.6) (el + C) k Fl e =6 u ,u£U, k £ N, 

in order to determine functions Fu . 

2. Let 5 W . V be the Kronecker delta. One has to solve \U\ linear system of the 
size \U\ x \U\ to determine coefficients a™ 

(2.7) *w,7 = £ «W,*(7). w^£U. 

3. It gives the following representation of the corresponding Lagrangian spline 

(2-8) L£ £ = £cCi^ £! w£U. 

ueu 

4. Every spline £ y(U, s, e) which takes prescribed values y = {y w }, w £ U, 
can be written explicitly as 

yj' = £ a,'-: ■ 

3. CUBATURE FORMULAS ASSOCIATED WITH VARIATIONAL INTERPOLATING 

SPLINES 

Theorem 3.1. In the same notations as above for every subset of vertices U = {u} 
and every k £ N, e > 0, there exists a set of weights 6 U — 8 U (U, k, s) £ R, u £ U, 
such that for every spline Y k £ that takes values Y k ' (u) = y u , u £ U, the following 
exact formula holds 

(3.i) £^» = £y<A 
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Proof. According to the above if is the Lagrangian spline , i. e. 



L k'e( V ) = 5 « 



for u G U, v G V, then the following representation holds 

(3.2) nT/(«) = E/( u K>)> veV - 



ueu 



By introducing scalars 



e u = e u {u, k,e) = Y, <'» 



□ 



For a subset S CV (finite or infinite) the notation L2(S) will denote the space 
of all functions from L2(G) with support in S: 

L 2 {S) = {ip G i 2 (G), = 0, v g 

Definition 1. We say that a set of vertices 5 C V is a A-set if for any ip g Li2(S) 
it admits a Poincare inequality with a constant A = A(5) > 

(3.3) M < A||ZV||, y> g L 2 (5). 

The infimum of all A > for which 5 is a A-set will be called the Poincare constant 
of the set 5 and denoted by A(S). 

The following Theorem gives a cubature rule that allows to compute the integral 
J2vev f( v ) ^ usm g on ly values of / on a smaller set U. 

Theorem 3.2. For every set of vertices U C V for which S — V \ U is a A-set 
and for any e>0,k = 2 l ,lGN, there exist weights 9 U = 9 U (U, k, e) g R such that 
for every function f g L^iG), 



(3.4) 



E /(») - E 



vEV 



ueu 



<2 v / \S\A k (eI + C) K f 



Proof. If / G Li2(G) and Y" fc ^ is a variational spline which interpolates / on a set 
U = V \ S then 



(3.5) 



E /(«) - E ifi'w 

Since 5 1 is a A- set we have 



< 



El/H-^w 



< 



f-K 



UJ 



< A 



II* (/-*£') 



(3-6) 

For any g G L^{G) the following inequality holds true 
(3-7) \\£g\\ < \\(sl + £)g\\. 

Thus one obtains the inequality 



(3.8) 



uj 

k.e 



< A 



{el + C) (/ 



-Y. 



uj 



k.e 



The following lemma holds true [3]. 
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Lemma 3.3. If A is a bounded self-adjoint positive definite operator in a Hilbert 
space H and for an tp £ H and a positive a > the following inequality holds true 

|M| <a\\A<p\\, 

then for the same ip £ H , and all k = 2 l , I — 0, 1, 2, ... the following inequality holds 

M<a k \\A k p\\. 

XJ f 

We apply this lemma with A = el + C, a = A and ip = f — Y k ' J ■ It gives the 
inequality 



(3.9) 



u,f 



< A k 



(eI + C) k [f-Y^ 



rU,} 



for all k — 2 l , / = 0, 1, 2, ... Using the minimization property of Y^'J we obtain 



f-y k , s 

Together with (|3.5j) it gives 



< 2A fc 



{el + Cf f 



(3.io) E/w-En^^ 

By applying the Theorem 13.11 we finish the proof. 



,k = 2\l G N. 



< 2 v / |^|A fc ||(e7 + /:) fe /||, = 2 ! ,i e N. 



□ 



The Bernstein inequality (|2.2[) and the last Theorem imply the following result. 

Corollary 3.1. For every set of vertices U C V for which S — V \ U is a A-set 
and for any e > 0, k = 2 l , I £ N, there exist weights 9 U — 9 U (U, k, e) £ R such that 
for every function f £ (£) , the following inequality holds 



(3.11) 



_ /(«) - E f^ 6 - 

vev ueu 



<2 7 fc /^ ll/H, 



w/iere 7 = A(w + e), k = 2 l , I £ N. 

If in addition the following condition holds 

1 

< oj < — - e 
A 

and / € E w (£) then this Corollary imply the following Theorem. 

Theorem 3.4. If U is a subset of vertices for which S = V \ U is a A-set then for 
any < e < 1/A, k = 2 l , I £ N, there exist weights 6 U = 8 U (U, k, e) £ R such that 
for every function f £ E U (C), where 



< oj < i - e, 



f/ie following relation holds 



(3.12) 

w/ien fc = 2' — > 00. 



E - E 



«6(7 
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Example 1. 

Consider the cycle graph Ciooo of 1000 vertices. The Laplace operator C as 
it defined in (I2.1[) has one thousand eigenvalues which are given by the formula 
A fc = 2-2cos^, fe = 0,l,...,999 (see [T]). 

It is easy to verify that every single vertex in Ciooo is a A = yg-set. It is also 
easy to understand that if closures of two vertices do not intersect i. e. 

(vj U dvj) n (vi U dvi) = 0, Vj, Vi G Ciooo, 
then their union Vj U Vi is also a A = -^-set. It implies, that one can remove from 
Ciooo every third vertex and on the remaining set of 670 the formula (|3.12[) will be 
true for the span of about 290 first eigenfunctions of C. 

Example 2. 

One can show [3] that if S = {v\, «2> — > Vn} consists of |5| < 998 successive 
vertices of the graph Ciooo then it is a A-set with 

A = -(sin— ^ — ) 

2\ 2\S\ + 2J 

It implies for example that on a set of 100 uniformly distributed vertices of Cioo 
the formula p. 121) will be true for every function in the span of about 40 first 
eigenfunctions of C. 

It is worth to note that the above formulas are optimal in the sense it is described 
bellow. 

Definition 2. For the given U C V, f G L 2 {G), k G N, e > 0, K > 0, the notation 
Q(U, /, k, e, K) will be used for a set of all functions h in L2(G) such that 

1) h(u) = f(u),ue u, 
and 

2) \\(eI + C) k h\\ < K. 

It is easy to verify that every set Q(U, /, fc, e, K) is convex, bounded, and closed. 
It implies that the set of all integrals of functions in Q(U, /, k,e,K) is an interval 
i. e. 

(3.13) [a,b] = lj2 h ^ : heQ(UJ,k,e,K)\ 

[v£V J 

The optimality result is the following. 

Theorem 3.5. For every set of vertices U C V and for any e > 0, k — 2 l , I G N, if 

(U, k, e) G M are the same weights that appeared in the previous statements, 
then for any g G Q(U, /, fc, e, K ) 

(3-14) J29(u)9 u = ^, 

where [a, b] is defined in (3.13\) . 



s 
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Proof. We are going to show that for a given function / the interpolating spline 

■ k.e 



Y^'J is the center of the convex, closed and bounded set Q(U, f, k,e, K) for any 



K > 



(eI + C) k Y 



k v UJ 
k.e 



. In other words it is sufficient to show that if 



Y, 



uj 



k.e 



h g Q(U,f,k,e,K) 



u f 

for some function h then the function Y k ' —h also belongs to the same intersection. 
Indeed, since h is zero on the set U then according to (|2.3j) one has 



(el + C) k Y^/, (el + C) k h) =((el+ C) 2k Y^/,h) = 0. 



But then 



(eI + C) k <yff + h) 



(eI + C) k (Y^ f -h 



In other words, 



(el + £) k ( 



Y 



uj 



< K 



and because YF'J + h and Y^'J — h take the same values on U the function Y h 
belongs to Q(U, f, k, e, K). From here the Theorem follows. 



rUJ 



UJ 

k,e 



□ 



Corollary 3.2. Fix a function f G Li2(G) and a set of vertices U C V for which 
S = V \ U is a A-set. Then for any e > 0, k = 2 l , I G N, for the same set of 
weights 9 U — 9 U (U, k, e) G ffi. that appeared in the previous statements the following 
inequalities hold for every function g G Q(U, /, k, e, K), 



(3.15) - E 

vev ueu 

Proof. Since / and g coincide on U from p. 51) and 



< V|S|A fc dmTO<3([/,/,fc,e,iO. 



we obtain the inequality 



(3.16) 



By the Theorem 13.51 the following inequality holds 



<v^A fc (el + Cflf-Y^ 



1 



for any g G Q(U, /, fc, e, if). The last two inequalities imply the Corollary. 



< —diamQ(U, /, k, e, K) 



□ 



4. CUBATURE FORMULAS FOR BANDLIMITED FUNCTIONS 

We introduce another set of cubature formulas which are exact on some sets of 
bandlimited functions. 

Theorem 4.1. If U is a subset of vertices for which S = V \ U is a A-set then 
there exist weights a u — a u (U) G M, u G U, such that for every function f G E U (C), 
where 

< uj < — , 
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the following exact formula holds 

(4.i) = u = v\s, 

vev ueu 

Proof. First, we show that the set U is a uniqueness set for the space E w (£), i. 
e. for any two functions from E u (£) the fact that they coincide on U implies that 
they coincide on V. 

If /, 5 € E w (£) then / — g G E w (£) and according to the inequality (|2.2j) the 
following holds true 

(4-2) \\C(f-g)\\<io\\f-g\\. 

If / and g coincide on U = V\S then f — g belongs to L 2 (S) and since S is a A-set 
then we will have 

||/ -. 9 || < A||£(/- 5 )||, f-geL 2 (S). 

Thus, if / — g is not zero and ui < 1/A we have the following inequalities 

(4-3) ||/ - ,g|| < A||£(/ - ,9)|| < Aw\\f - g\\ < \\f ~ g\\, 

which contradict to the assumption that / — g is not identical zero. Thus, the set 
U is a uniqueness set for the space E w (£). 

It implies that there exists a constant C = C(U, ui) for which the following 
Plancherel-Polya inequalities hold true 

1/2 / \ 1/2 

(4-4) ' " 



Ei/^i 2 ) < ii/ii <cfEi/(")i 2 ) 

sueu I \u£U I 



for all / £ E w (£). Indeed, the functional 

/ \ 1/2 

111/111= E I/hi 



\u£U / 

defines another norm on E w (£) because the condition |||/||| = 0, / G E u (£), implies 
that / is identical zero on entire graph. Since in finite-dimensional situation any 
two norms are equivalent we obtain existence of a constant C for which (|4.4[) holds 
true. 

Let 8 V G L 2 (G) be a Dirac measure supported at a vertex v s V. The notation 
will be used for a function which is orthogonal projection of the function 

1 . 



on the subspace E^(£). If ipo, tpi, (pjrw) are orthonormal eigenfunctions of £ 
which constitute an orthonormal basis in E w (£) then the explicit formula for § v is 

(4.5) i = ^ft'WW- 

3'=0 

In these notations the Plancherel-Polya inequalities (|4.4j) can be written in the form 

(4.6) Ei^^i'-ii/ii^^Ei^^)! 2 ' 

ueu ueu 

where /, d u £ E w (£) and (/,i?«) is the inner product in L 2 (G). These inequalities 
mean that if U is a uniqueness set for the subspace E w (£) then the functions 
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{$u}u£U form a frame in the subspace E w (£) and the tightness of this frame is 
1/C 2 . This fact implies that there exists a frame of functions {& u }ueu in the space 
E w (£) such that the following reconstruction formula holds true for all / <E E w (£) 



Unfortunately this approach does not give any information about constant C in 
(I4.6|) and it make realization of the Theorem l4 . 1 [ problematic . We are going to utilize 
another approach to the Plancherel-Polya-type inequality which was developed in 
our paper [7] and which produces explicit constant. 

For any U which is a subset of vertices of G we introduce the following operator 



(4.8) d°(U) = U, cl(U) = SUdU, d rn (U) = d (cJ m_1 (I0) . m e N, U C V. 



We will use the following notion of the relative degree. 

Definition 3. For a vertex v £ d m (U) we introduce the relative degree d m (v) as 
the number of vertices in the boundary d (cl m (U)) which are adjacent to v. 



Definition 4. For a vertex v € d (d m (U)) we introduce the quantity k m (v) as the 
number of vertices in the set d m (U) which are adjacent to v: 



(4.7) 




By setting a, 



J2vev one obtains (|4TJ) . 



□ 



d m (v) =card{w Ed{d m (U)) : w ~ «} . 
For any S C F we introduce the following notation 

-Dm = D m (U) = sup d m (v). 

t)6cZ m ([/) 



k m (v) = card{w e d m (U) : w ~ u} . 



For any U G V we introduce the following notation 



K m =K m (U)= inf fcm(v). 

t)G9(ci m (C/)) 



For a given set U CV and a fixed n S N consider a sequence of closures 



U,cl(U),...,cl n (U),n € N. 



Theorem 4.2. in i/ie sa?7ie notations as above if the condition 
(4.9) cZ"(t/) = F 

and i/ie inequality 



(4.10) 




/ioW, t/ien /or any / £ -EL (£) t/ie next inequality takes place 




(4.11) 




1/2 
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where 

Theorem 4.3. If U is a subset of vertices for which condition J^.ffp is satisfied 
then there exists a set of weights ji u 6 R, u € £/, smc/i i/iat /or any / £ E^C), 
where uj satisfies {4-.10fy the following exact formula holds 

(4.i3) a«) = E 

uev net/ 

Proof. The previous Theorem shows that J7 is a uniqueness set for the space E u (£), 
which means that every / in E w (£) is uniquely determined by its values on U. This 
fact implies existence of a frame {</>«}, u £ U, in E w (£) such that 

, / \ 1/2 

£k/,<mi 2 < ii/ii < 



1-7 



1 -7 v 

' \i=0 v 7 / \«€J7 



1/2 / v 1/2 



2 



where 7 as in (|4.12j) . 

Now, by using a frame which is dual to the frame {4>u} and by setting 

M« — S-ugv we obtain the Theorem. 

□ 

All the results in this section are only useful if it is possible to calculate dual 
frame explicitly. Since it is not a simple task we are going to offer a certain way 
around which is based on the so-called frame algorithm [2], Ch. 5. 

Let {ej} be a frame in a Hilbert H space with frames bounds A, B, i. e. 

A||/|| 2 <^|</, ei )| 2 <B||/|| 2 , fen. 

Given a relaxation parameter < A < -g, set 8 = max{|l — XA\, |1 — XB\} < 1. 
Let fo = and define recursively 

(4.15) /„ = /„_! + XS(f - /„_!), 

where S is the frame operator which is defined on H by the formula 

3 

In particular, f\ — XSf = XJ2j (f> e j) e j- Then limn-^ /„ = / with a geometric 
rate of convergence, that is, 

(4-16) ||/-/n]]<<H|/||. 

Note, that for the choice A = A ^_ B the convergence factor is 
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Let us go back to our situation in which we consider a subset of vertices U C V. 
Just for the sake of simplicity we will consider the case when 

(4.18) U = UUdU = V. 

It can be useful to remind that in this case the relative degree do(v) is the number 
of vertices in the boundary dU which are adjacent to a v € U: 

do(v) = card{w G dU : w ~ v} 

and 

D = D {U) = sup d (v). 

For a vertex v G dU, the quantity ko (v) is the number of vertices in the set U which 
are adjacent to v: 

ko(v) — card{w G U : w ~ v} 

and 

X = K (U) = inf fc («). 
We consider the Hilbert space H is E w (£), where 

(4.19) u < 

The set 17 is a uniqueness set for the space E w (£) and the corresponding frame in 
it is {4> u }, u G U. The frame operator S takes the following form 

(4.20) Sf=Y,(fAu)<l>u, /eE u (4 

and 

(f,ct> u )=f(u), f G E w (£), u E U. 
Thus the recurrence sequence (|4.15l) takes the form /o = 0, and 

(4-21) /„ = /„_a + A ^(/ - f n -i)(u)d> u , 

ueu 

We are ready to state the following fact which provides a reconstruction method 
for functions in E w (£) from their values on {7. 

Theorem 4.4. Under the assumptions . _? /or / G E u (£) the following 

inequality holds for all natural n 

(4-22) ||/- /„|| < S»\\f\\, 

where the convergence factor is 



8 



D (U) + K (U)- 
Proof. The frame inequality (|4.2I) takes the form 

where / G E w (£), /(u) = (/, U ) , u <E U, and 



7 = 2y/u/K (U) < 1. 
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It shows that the frame bounds of the frame {4> u } arc 

^2 K (U) 



A = (i- 1 y 



B = (l-7) • 



2D (U) + K (UY 

For /€E U (£) and /„ defined in (|4. 15)) one has the inequality (14. 16)) and according 
to (|4.17p the convergence factor is 

D (U) 



(4.24) 

The Theorem is proved. 
Since 



8 



~ D (U) + K (Uy 
D a (U) 



□ 



Do(U)+K (U) l + 

it becomes clear that the best convergence occurs for such sets U which maximize 
the ratio of the number K (U) of "in going" edges to the number D (U) of "out 
going" edges. 
By setting 

v u = E <Pu(v), u elf, 

v£V 

we obtain the recurrence sequence for corresponding "integrals" 



E /o(«) = o> 
E/i(«) = aE/( u K> 

v£V u£U 

E = E + A E(/ - u-x){u)v u . 



vev vev ueu 

Thus, the "integral" ^2 v£V f n (v) for every f n is expressed in terms of weights 
{^ u }, u <E U, and values of / on the uniqueness set U. 

We are ready to state the following fact which provides an approximate cubature 
formula for functions in E w (£) . 

Theorem 4.5. Under the assumptions \4- 18\ ), for f G E W (C) andJ2 v evf n ( v ) 

defined in (|4]) we have 



E/( u )-E/«w 



vev 



vev 



< 8 n 



where the convergence factor is 



5 = 



D Q (U) 



D (U)+K (U)' 



The proof is obvious: 



<El/( V )-/«( W l^ VW\\\f-fn\\<S 



vev 



E/h-Ez-w 

vev vev 
Example 3. 

If one will consider the cycle graph Ciooo and will remove from it every third 
vertex then the set U will contain about 670 remaining vertices and the space E w (£) 
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will have dimension 290 (see Example 1). In this case Dq(U) = 1, Kq(U) = 2 and 
the convergence factor S = 1/3. 
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